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FOREWORD 


The  latest  version  of  the  NSWCDD  Aeroprediction  code  (AP95)  produces  linear 
aerodynamic  loads  that  are  based  on  theoretical  methods  which  compute  pressures  all  along  the 
body  and  lifting  surfaces.  These  linear  loads  can  be  integrated  to  obtain  total  shear  and  bending 
moments  for  use  in  structural  design  and  analysis.  The  AP95  also  computes  the  total  nonlinear 
aerodynamic  loads  but  not  in  pressure  distribution  form.  Thus,  the  nonlinear  term  of  the 
aerodynamic  load  (which  is  quite  large  at  high  angle  of  attack)  is  not  included  in  the  shear  and 
bending  moment  calculations.  As  a  result,  the  AP95  is  not  as  useful  to  the  structural  design 
engineer  as  it  could  be.  To  remedy  this  situation,  new  technology  was  developed  to  distribute 
these  nonlinear  loads  on  the  various  missile  or  projectile  components.  This  new  technology  is 
discussed  in  this  report. 
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1.0  INTRODUCTION 


The  Naval  Surface  Warfare  Center,  Dahlgren  Division  (NSWCDD)  Aeroprediction  Code 
(AP95)^  allows  the  computation  of  aerodynamics  over  an  extensive  range  of  flight  conditions 
which  includes  the  operating  environments  encountered  by  most  tactical  missiles.  This  range 
covers  angles  of  attack  (AOA)  up  to  90  deg,  Mach  numbers  from  0  to  15,  and  control  deflections 
from  -30  to  +30  deg.  Originally,  it  was  hmited  to  the  roll  position  of  O  =  0  deg  (fins  in  the  plus 
(+)  position),  but  Reference  2  recently  extended  the  code  to  the  d)  =  45  deg  (fins  in  the  cross  (x) 
position)  roll  orientation. 

The  latest  public  release  of  the  AP95  code  incorporates  the  capability  to  model  nonlinear 
aerodynamic  loads.  These  nonlinear  loads  are  estimated  in  an  approximate  analytical  or 
semiempirical  sense  by  separating  each  aerodynamic  force  component  into  a  linear  and  nonlinear 
term.  As  in  previous  versions  of  the  code,  the  linear  term  is  modeled  by  the  application  of  either 
linear,  slender  body,  or  second-order  theory.  The  resulting  pressure  distributions  can  be 
integrated  to  produce  aerodynamic  loading  information  for  use  in  structural  design  and  analysis. 
The  nonlinear  contribution,  however,  is  modeled  by  the  direct  application  of  wind  tunnel  data 
bases  and  is  defined  primarily  in  terms  of  total  force  and  moment  measurements.  Some 
distribution  information  is  provided  for  the  nonlinear  component  of  the  body  loads,  but  none  is 
currently  available  for  the  lifting  surface  loads  or  for  their  interference  effects  on  the  body.  This 
situation  is  not  restrictive  in  terms  of  aerodynamic  and  performance  analysis.  Structural 
engineers,  however,  need  to  know  not  only  the  magnitude  of  the  aerodynamic  forces,  but  also 
how  they  are  distributed  over  the  surface  of  a  missile  if  they  are  to  determine  the  shear  and 
bending  moments  to  which  its  components  will  be  subjected.  For  this  reason,  the  AP95  code  is 
being  modified  to  allow  the  prediction  of  the  distribution  of  the  nonlinear,  as  well  as  the  linear, 
aerodynamic  loads  over  both  the  body  and  control  surfaces.  All  interference  effects  wiU  be 
included  in  the  analysis  in  addition  to  the  individual  component  nonlinear  aerodynamics.  It  is 
believed  that  the  next  version  of  the  Aeroprediction  code  (AP98)  will  be  the  first  code  of  its  kind 
to  offer  this  capability. 


2.0  ANALYSIS 


The  total  normal  force  for  a  wing-body-tail  configuration  can  be  defined  by: 


C)f=Cfj  +AC;y  +AC^  +C 


’£(7)  ^^T(V) 


(1) 
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The  &st  term  in  Equation  (1)  is  the  normal  force  coefficient  of  the  body  alone.  The  second  term  is 
the  normal  force  coefficient  of  the  wing  in  the  presence  of  the  body  and  the  third  represents  the 
additional  component  of  normal  force  on  the  body  due  to  the  presence  of  the  wing.  The  fourth 
term  is  the  normal  force  coefficient  of  the  tail  in  the  presence  of  the  body  and  the  fifth  represents 
the  incremental  normal  force  on  the  body  as  a  result  of  the  tail.  The  final  term  is  the  normal  force 
on  the  tail,  usually  negative,  caused  by  the  vortices  shed  from  the  wing. 

Each  of  the  terms  in  Equation  (1)  has  a  linear  and  a  nonlinear  component.  Either  linear 
theory  or  slender  body  theory  is  used  to  determine  the  linear  portion  of  each  term  and  the 
nonlinear  contributions  are  computed  directly  from  wind  tunnel  data  bases  with  engineering 
judgment  and  extrapolation  being  used  where  data  base  information  is  limited.  Overall  average 
accuracy  is  maintained  at  ±10  percent  for  the  total  normal  and  axial  force  coefficients  and 
±4  percent  of  total  body  length  for  the  center  of  pressure. 

In  the  following  sections,  each  component  of  Equation  (1)  wiU  be  considered  individually 
and  the  methodology  which  has  been  developed  to  distribute  its  nonlinear  elements  over  the  body 
surface  will  be  presented.  First,  the  d>  =  0  deg  roll  position  wiU  be  addressed  and  then  any 
changes  or  additions  which  are  necessary  to  extend  the  methodology  to  the  <l>  =  45  deg  roll 
position  will  be  presented. 


2. 1  BODY  LOADS  AT  THE  <1)  =  0  DEG  ROLL  POSITION 

Four  of  the  terms  in  Equation  (1)  contribute  to  the  body  loading.  These  are  , 
AC„  ,  ACm  ,  and  ■  Only  the  first  three  of  these  will  be  discussed  in  this  section.  The 

wing-tail  interference  will  be  reserved  for  the  section  on  lifting  surface  loads  since  it  is  connected 
directly  to  them  and  is  shared  with  the  tail  surfaces. 


2.1.1  Body  Alone  Loads 


The  body  alone  component,  ,  is  composed  of  a  linear  and  a  nonlinear  load.  The  linear 

load  is  computed  analytically  above  M  =  1.2  by  either  Van  Dyke’s  hybrid  theory^  or  second  order 
shock-expansion  theory.'*  Pressure  coefficients  are  computed  as  a  function  of  position  along  and 
around  the  body  and  then  integrated  along  the  body  to  obtain  local  linear  normal  force.  That  is 


)i. 


da 


dx 


-2 


Tvr: 


jCp{x,  <l>)(cos  ^)rd^ 


(2) 


The  nonlinear  body  normal  force  coefficient  contribution  is 

^REF 


(3) 
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If  we  assume  that  the  nonlinear  load  is  distributed  uniformly  along  the  body  surface,  then  the  local 
nonlinear  normal  force  contribution  at  any  point  x  is 


7]Co  sin^a 


•(x)  . 


^REF 


(4) 


The  total  local  normal  force  coefficient  at  any  point  along  the  body  surface  is  the  sum  of 
Equations  (2)  and  (4). 


Comparison  of  Equation  (5)  with  results  from  Computational  Fluid  Dynamics  (CFD) 
computations  indicated  that  the  predicted  body  load  in  the  nose  region  was  too  high  at  low  AOA. 
It  is  beheved  that  this  discrepancy  occurs  because,  at  low  AOA,  the  point  at  which  the  body 
vortices  begin  to  separate  is  toward  the  afterbody  of  the  missile.  Then,  as  AOA  is  increased,  the 
separation  point  moves  forward  and  the  full  effect  of  the  body  shed  vortices  is  apparent.  In  order 
to  account  for  this  effect,  the  nonlinear  load  contribution  is  adjusted  as  follows: 

X  ^  LfioSE  ’ 


NL,ADJ 


NL 


^NOSE  X:  <  2D 


'HOSE 


^”b)nL,ADJ  r 


X-  L. 


'NOSE 


\  j 


(6) 


(7) 


X  ^  Linose  2Dj  : 


^NL,ADJ  )nL  [(^”*  )nL  ]iV05£  irflAIV 

On  the  nose  itself,  as  shown  in  Equation  (6),  the  nonlinear  load  is  reduced  at  AOA  below 
a  =  30  deg  by  a  factor  which  varies  linearly  with  AOA  from  1  at  a  =  30  deg  to  0  at  a  =  0  deg. 
No  adjustment  is  made  at  AOA  above  30  deg.  There  is  a  transition  region,  which  extends  for  a 
distance  of  two  calibers  from  the  end  of  the  nose,  within  which  the  nonlinear  load  component  is 
allowed  to  recover  gradually  from  its  reduced  value  at  the  end  of  the  nose  to  its  normally 
computed  level  at  the  current  location,  as  shown  in  Equation  (7).  Since  the  total  normal  force 
computed  by  AP95  agrees  well  with  experimental  results,  the  portion  of  the  nonlinear  loads  which 
are  removed  from  the  nose  and  transition  regions  must  be  added  back  into  the  loads  on  the 
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remainder  of  the  body.  Thus,  Equation  (8)  shows  that,  for  points  on  the  body  more  than  two 
calibers  downstream  of  the  nose,  the  normally  computed  local  nonlinear  loads  are  increased  by 

two  terms  I  (c„  ]  ]  and  ffc„  ]  1  which  represent  the  total  load  reductions  in  the  nose 

and  transition  regions,  respectively,  spread  evenly  over  the  total  remaining  body  length.  Figure  1 
shows  the  effect  of  this  change  on  the  body  loads  for  a  typical  case,  using  the  SEASP ARROW 
air-to-air  missile  as  an  example.  As  can  be  seen,  relative  to  the  CFD  results  which  are  shown  for 
comparison,  there  is  some  improvement  in  the  nose  region  and  little  change  elsewhere. 


2.1.2  Bodv-Wing  and  Body-Tail  Interference  Loads 

For  the  discussion  to  follow,  forward  lifting  surfaces  will  be  referred  to  as  wings  and  aft 
lifting  surfaces  will  be  referred  to  as  the  tail.  In  the  vicinity  of  any  lifting  surface,  the  body  is 
subjected  to  increased  loads  caused  by  interference  effects.  These  contributions,  represented  in 
Equation  (1)  by  for  the  wing  and  by  for  the  tail,  are  also  computed  as  separate 

linear  and  nonlinear  components  with  afterbody  effects  being  considered  for  the  rear  lifting 
surfaces.  For  simplicity,  it  was  first  assumed  that  these  interference  loads  were  distributed 
uniformly  along  the  root  chord  of  the  wing  or  tail,  thus  affecting  the  body  loading  only  along  the 
attachment  line  from  the  leading  edge  to  the  trailing  edge.  These  constant  local  loading  values 
were  defined  as 


and 


Ac„ 


AC. 


(cX 


Ac„._.  =■ 


AC^ 


"B(T) 


(9) 

(10) 


Figure  2(a)  shows  the  region  of  the  body  affected  by  these  loads  and  Figure  2(b)  shows  the 
distribution  of  their  magnitude.  Comparisons  with  CFD  data  indicated  that  this  approach 
produced  interference  effects  which  were  too  localized,  with  excessively  high  values  along  the 
root  chords  and  none  of  the  carryover  load  behind  the  lifting  surfaces  indicated  by  the  data. 
Therefore,  the  region  of  influence  of  the  lifting  surfaces  was  changed  to  coincide  with  the  area  on 
the  body  surface  bounded  by  the  Mach  lines  coming  from  their  leading  and  trailing  edges,  as 
shown  in  Figure  3(a).  The  new  distribution  of  the  load  magnitude  on  the  body  is  shown  in  Figure 
3(b).  It  is  consistent  with  the  method  currently  used  for  moment  and  center  of  pressure 
calculations  in  the  AP95.  Note  that  the  magnitude  increases  linearly  from  the  leading  edge  to  the 
point  where  the  Mach  lines  from  opposite  lifting  surfaces  meet  in  the  center  of  the  body.  Since 
the  angle  of  the  Mach  fines,  q,  is  given  by 


fi  =  sin 


(11) 
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the  streamwise  location  where  this  intersection  occurs  is 


^2 


(12) 


where 


P  =  ylMl-l  .  (13) 

The  local  Mach  number,  Ml,  is  approximated  by  the  freestream  Mach  number.  For  the  region  on 
the  body  enclosed  by  both  Mach  lines,  the  entire  body  surface  is  affected  by  the  interference  and 
the  load  value  is  held  constant  at  the  maximum  level.  At  the  streamwise  location  given  by 

x,=x^+^  +  C,  (14) 


the  Mach  lines  from  the  trailing  edges  intersect  and  a  portion  of  the  body  begins  to  emerge  from 
the  interference  region  until,  at  the  point  defined  by 

X4=x^+I3Db+C,  (15) 


aU  Mach  lines  have  moved  off  the  body  and  the  interference  loads  are  reduced  to  zero.  Between 
X3  and  X4,  the  interference  load  is  allowed  to  decay  linearly  from  its  maximum  value  to  zero. 


Since  the  total  value  of  the  interference  normal  force  contribution  must  remain  the  same  as 
for  the  uniform  distribution  originally  assumed,  it  is  necessary  to  adjust  the  maximum  value  of  the 
local  interference  load  relative  to  the  constant  value  of  Equations  (9)  and  (10)  and  Figure  2(b). 
Since  both  the  wing  and  the  tail  adjustments  are  handled  in  an  identical  manner,  the  W  and  T 
subscripts  wiU  be  dropped  and  only  one  set  of  relations  will  be  shown.  The  maximum  values  will 
become 


2(q)  ^ 

^2{C,)  +  PD,j 


(16) 


Thus,  a  trapezoidal  load  pattern  is  produced,  defined  by 


Xle  ^x<x^: 


(17) 
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X2^x<Xj: 


(18) 


X2<x<x^: 


(19) 


In  situations  where  the  wing  or  tail  is  located  near  the  aft  end  of  the  body,  it  is  possible  for  the 
Mach  lines  to  extend  beyond  the  body  and  some  of  the  interference  loading  will  be  lost  unless 
further  adjustments  are  made.  For  these  cases,  the  cross-sectional  area  of  the  body  enclosed  by 
the  Mach  lines  and  the  aft  end  of  the  body  is  determined  and  compared  to  the  total  cross- 
sectional  area  that  would  be  enclosed  by  the  Mach  lines  if  the  body  were  not  cut  off.  The 
maximum  interference  load  values  are  then  adjusted  to  compensate  for  the  lost  area  and  maintain 
the  correct  total  of  interference  normal  force  as  computed  by  the  AP95. 


2.2  LIFTING  SURFACE  LOADS  AT  THE  0  =  0  DEG  ROLL  POSITION 

The  total  contributions  to  the  normal  force  of  the  wing  and  tail  are  represented  in 
Equation  (1)  by  and  .  respectively.  Using  as  an  example,  these  two  terms 

each  have  the  following  form; 

~  [•^w'(5)  sin  a  +  sin5j^C;^^  (20) 

Nonlinearities  are  introduced  into  by  defining  both  a  linear  and  a  nonlinear  contribution  for 

the  normal  force  of  the  isolated  wing,  represented  by  ,  and  for  each  of  the  interference 

factors  and  which  are  introduced  to  account  for  the  effects  of  the  presence  of  a  body 

on  the  normal  force  of  the  lifting  surface.  The  method  which  has  been  developed  to  distribute  the 
nonlinear  contributions  to  the  wing-alone  normal  force  and  the  interference  factors  wiU  be 
presented  in  the  next  section. 

The  final  term  in  Equation  (1)  ,  ,  is  also  associated  with  the  lifting  surfaces.  It 

represents  the  impact  on  the  tail  of  vortices  shed  from  the  wing  when  both  sets  of  lifting  surfaces 
are  present.  A  portion  of  this  downwash  effect  should  be  allocated  to  the  body,  since  it  also 
experiences  the  effect  of  the  shed  vortices.  Techniques  to  accomplish  both  the  body/tad  division 
of  this  load  contribution  as  weU  as  its  spanwise  distribution  along  the  tad  surfaces  wdl  be 
presented  below. 
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2.2.1  Nonlinear  Wing  or  Tail  Alone  and  Interference  Factor  Load  Distributions 

The  linear  part  of  the  local  normal  force  coefficient  on  a  wing  or  tail  surface  at  a  spanwise 
location  y  is 


l^’’{y)]i,=-^j^Cp{x,y)dx  .  (21) 

ACp(x,y)  is  the  difference  between  the  pressure  on  the  upper  and  lower  surfaces  at  a  given  (x,y) 
location.  This  pressure  differential  is  computed  from  either  lifting  surface  theory  for  subsonic 
flows  or  3-D  thin  wing  theory  for  supersonic  flows.  Equation  (21)  gives  the  spanwise  distribution 
of  the  linear  contribution  to  the  normal  force  of  an  isolated  lifting  surface.  The  total  linear 
contribution  to  normal  force  can  be  found  by  integrating  Equation  (21)  in  the  spanwise  direction. 

The  nonlinear  part  of  the  wing  alone  lift  is  determined  using  a  fourth  order  analytical 
method^  in  conjunction  with  wing  alone  data  bases.^’^’^  Only  the  total  values  of  lift  and  moment 
are  available  from  the  data  bases  with  no  surface  distribution  information  included.  This 
distribution  wfll  be  determined  in  conjunction  with  that  of  the  interference  factors. 

From  the  combined  linear  and  nonlinear  wing-alone  normal  force,  is  determined 

based  on  a  secant  slope  of  the  total  angle  of  attack,  |  a-i-5 1 ,  experienced  by  the  lifting  surface. 
This  slope  is  then  multiplied  by  sin  a  and  sin  6  terms  which  contain  the  interference  factors  K  and 
k,  respectively,  to  account  for  the  interference  effects  of  the  body  on  the  wing.  (All  subscripts 
have  been  dropped  from  K  and  k  for  simplicity.)  The  linear  contributions  to  K  and  k  are  found 
from  slender  body  theory.  These  linear  values  are  then  adjusted  for  nonlinear*^’^’^  effects  based  on 
information  from  wind  tunnel  data  bases.  As  for  the  wing-alone  case,  these  data  bases  provide 
only  total  force  and  moment  measurements  and,  therefore,  a  method  must  be  developed  to 
distribute  the  nonlinear  component  of  K  and  k  along  the  span. 

As  shown  above  in  Equation  (21),  the  AP95  code  already  computes  the  spanwise 
distribution  of  the  linear  normal  force  of  an  isolated  lifting  surface.  It  will  be  assumed  that  the 
nonlinear  contributions  to  the  wmg-alone  normal  force  and  to  the  interference  effects,  represented 
by  the  K  and  k  factors,  wiU  be  distributed  in  the  same  manner.  Thus,  at  any  y  location,  the  local 
load  c„(y)  for  the  wing  and  tail  will  be 


(c„(y))H,  =([c„(y)]i)^ 


C. 


^W(S) 


\SPAN 


(22) 


and 
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c. 


X([^.W]J 

\SPAN 


W 


(23) 


2.2.2  Wing-Tail  Interference 

The  AP95  code  computes  only  the  total  value  of  the  wmg-tail  interference  term  ( in 

Equation  (1)).  In  this  instance,  however,  not  only  is  a  method  of  spanwise  distribution  on  the  tail 
needed,  but  it  is  also  apparent  that  a  portion  of  the  downwash  effects  should  be  carried  over  to 
the  body.  The  division  of  between  the  tail  and  body  will  be  considered  first.  The  method 

in  which  this  objective  is  accomplished  is  not  rigorous,  but  does  have  some  intuitive  foundation 
and  is  probably  as  well  as  can  be  done  with  the  limited  amount  of  CFD  data  available. 

Refer  to  Figure  4  for  the  geometry  of  this  model.  The  assumption  is  made,  from  slender 
body  theory,  that  the  vortices  from  the  forward  lifting  surfaces  are  shed  from  a  spanwise  location, 
yvORTj  of 


yvORT  ~ 


(24) 


and  travel  straight  back  parallel  to  the  velocity  vector  from  this  point.  A  relative  value  for  the 
effect  of  this  shed  vortex  on  the  tail  surface  at  any  spanwise  point  is  given  by 


VC{y)  = 


1 


^VORT  yDIST 


(25) 


where 


y  BIST  ~  lyVORT  H  (^^) 

and  rvoRT  is  a  “vortex  radius”  taken  to  be  0.04.  This  value  of  rvoRi  could  be  subject  to  adjustment 
in  the  future  as  additional  CFD  or  experimental  data  becomes  available.  An  inspection  of 
Equation  (25)  shows  that,  at  any  point  y,  it  consists  of  the  product  of  a  factor  which  is 
proportional  to  the  local  vortex  strength  through  the  inverse  square  distance  term  and  a  second 
factor,  represented  by  the  local  chord,  which  is  related  to  the  local  area  affected  by  the  vortex. 
The  use  of  an  absolute  value  for  yoisT  assures  that  the  relation  will  be  valid  regardless  of  where  the 
vortex  location  falls  on  the  tail  span,  even  if  it  lies  outboard  of  the  tip.  The  function  of  rvoRx  is  to 
prevent  singularities  if  this  distance  becomes  zero.  On  the  body,  a  similar  relation  is  used. 
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VC(y)  = 


^  1  Y 

VVORT  yOIST  > 


(27) 


Note  that  the  use  of  the  root  chord  for  the  area  factor  assumes  that  only  the  body  area  directly 
between  the  tail  surfaces  is  affected.  This  assumption  is  not  strictly  correct,  but  it  offers  a  starting 
point  until  further  refinements  can  be  made. 


Using  these  relations,  the  apportionment  of  the  wing-tail  interference  losses  can  be 
accomplished  by  summing  the  VC  terms  for  both  the  body  and  tail  and  then  taking  the  ratios  of 
the  stimmations  to  their  combined  total. 


BODY 


lVC(y)  ] 

_ BODY _ 


\BODY  TAIL  J 


(28) 


\c  1  =  c 

'lyci.y)  t 

TAIL 

L  \tail 

E^c(>')+E^^(>’) 

^  BODY  TAIL  y 

(29) 


It  will  be  assumed  that  the  portion  of  the  wing-tail  interference  load  which  is  allocated  to  the  body 
by  Equation  (28)  is  distributed  evenly  along  the  body  directly  between  the  tail  surfaces. 


(30) 


This  additional  contribution  is  added  to  the  total  local  body  load  at  streamwise  locations  between 
the  leading  and  trailing  edges  of  the  tail  surfaces. 


The  portion  of  the  wing-tail  interference  loads  assigned  to  the  tail  surfaces  by 
Equation  (29)  is  distributed  across  the  span  of  the  tail  by  assuming  that  its  magnitude  at  any  y 
location  is  proportional  to  the  VC(y)  value  (from  Equation  (25))  at  that  location. 


(^n  (y ))7’(V) 


(  ^ 

K  1 

L  "r(V)  JrA/i 

VC{y) 

JyCiy) 

\SPAN  y 

(31) 
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These  loads  are  added  to  those  computed  in  Equation  (23)  to  get  the  total  local  load  at  any 
spanwise  location. 


2.3  CHANGES  FOR  THE  O  =  45  DEG  ROLL  POSITION 

For  the  O  =  45  deg  roll  position,  several  changes  must  be  made  to  the  above  methodology 
to  produce  a  more  accurate  representation  of  the  surface  loads.  It  is  obvious  that  there  must  be 
some  way  of  accounting  for  the  different  aerodynamic  environments  experienced  by  the  windward 
and  leeward  plane  fins.  In  addition,  the  distribution  of  the  body  carryover  loads  induced  by  the 
lifting  surfaces  behave  in  a  somewhat  different  fashion  than  for  the  <D>  =  0  deg  case.  The  following 
sections  describe  the  modifications  that  were  made  to  the  methodology  described  above  to 
accommodate  these  differences. 


2.3.1  Bodv-Wing  and  Body-Tail  Interference  Load  Changes 

The  AP95  body  alone  normal  force  is  considered  to  be  independent  of  roll  position. 
However,  the  body-wing  and  body-tail  interference  effects  will  change  with  the  roll  configuration. 
The  O  =  45  deg  interference  factors  are  computed  automatically  when  this  option  is  selected  in 
the  input  file  and  are  simply  substituted  for  and  in  aU  computations.  The  values 

of  these  interference  factors  were  chosen  to  give  good  agreement  with  experiment  for  the  total 
loads  and  center  of  pressure.  Unfortunately,  when  compared  to  CFD  data,  it  appears  that  they  do 
not  produce  an  accurate  distribution  of  those  loads  along  the  body  or  on  the  leeward  plane  tail 
fins.  There  appear  to  be  three  physical  phenomena  which  account  for  these  discrepancies.  The 
first  of  these  will  be  discussed  now  and  the  others  in  the  following  section. 

On  the  surface  of  the  body  between  the  Ufting  surfaces,  there  is  a  region  of  high  dynamic 
pressure.  Immediately  behind  the  lifting  surfaces,  however,  if  there  is  an  afterbody  present,  there 
is  an  area  where  the  dynamic  pressure  is  much  lower.  The  current  AP95  methodology  at 
O  =  45  deg  averages  this  difference  out  to  get  the  right  total  normal  force  and  adds  a  center  of 
pressure  shift  which  accounts  for  the  inconsistencies  in  distribution.  If  the  distribution  of  the  body 
carryover  loads  is  to  be  modeled  correctly,  some  consideration  must  be  given  to  the  physics  of  the 
problem.  What  appears  to  be  happening,  based  on  the  experimental  data,‘  is  that,  in  the  region  of 
the  Ming  surfaces,  there  is  an  enhancement  to  the  body  load  at  both  the  <I)  =  0  deg  and  <I>  =  45 
deg  roll  positions.  At  ^>  =  0  deg,  this  enhancement  can  extend  several  body  diameters 
downstream  into  the  traditional  Mach  line  influence  region.  The  method  employed  to  model  this 
effect  was  discussed  in  Section  2.1.2.  At  fl>  =  45  deg,  the  body  loads  are  influenced  to  a  lesser 
extent  in  the  immediate  vicinity  of  the  Ming  surfaces,  but  are  still  much  higher  than  in  the  body 
alone  case.  Downstream  of  the  fins,  however,  the  high  dynamic  pressure  appears  to  decrease 
rapidly  and  produce  a  region  of  overexpansion  where  the  body  carryover  actually  becomes 
negative  and  total  loads  fall  below  those  of  the  body  alone.  Thus,  with  an  afterbody  present,  the 
body  carryover  load  is  lower  than  with  no  afterbody.  (Note  that  the  reverse  is  true  at  the 
0  =  0  deg  roll  orientation.) 
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Reference  2  models  this  effect  with  a  lower  minimum  value  of  Kb(w)  at  O  =  45  deg  versus 
<I>  =  0  deg  and  a  forward  shift  of  the  center  of  pressure.  This  approach  does  not,  however, 
produce  the  proper  distribution  of  the  interference  loads.  For  the  present  work,  the  following 
technique  was  used  to  accomplish  this  goal.  It  is  first  assumed  that,  for  the  purpose  of  computing 
the  body  carryover  load  distribution,  the  minimum  value  of  Kb(W)  is  the  same  at  0  =  45  deg  as  it 
is  at  =  0  deg.  This  value  of  [Kb(W)]min  is  distributed  in  the  vicinity  of  the  lifting  surfaces  for 
both  roll  orientations  as  described  in  Section  2.1.2.  However,  for  the  O  =  45  deg  roll  position,  an 
amount  equal  to 


is  subtracted  from  Kb(W)  in  the  region  downstream  of  position  X4  where  X4  is  defined  by  Equation 
(15).  Figure  5  shows  how  this  negative  load  is  distributed.  It  is  assumed  that 

\Ac„  ]  =-[Ac„  1  (33) 

L  "b  IMIN  L  ”b \mAX 

and  that  the  slope  from  X4  to  xs  is  the  same  as  the  slope  from  X3  to  X4.  The  slope  from  xg  to  x?  is 
the  negative  of  this  value.  The  positive  contribution  for  <I)  =  0  deg  has  already  been 

established  and  so  the  total  negative  contribution  can  be  found  as 

(34) 

where  the  W  subscript  has  been  dropped  to  avoid  repetition.  This  relation  and  those  to  follow 
win  be  the  same  for  both  wing  and  tail  surfaces.  It  can  be  seen  that  the  total  value  of  the  body 
carryover  term  is  not  changed.  The  other  defining  parameters  can  be  computed  as  follows: 


Xj  =  X4  +  (X4  -  X3  ) 

(35) 

(36) 

Ac;^^  =Ac;J^ 

(37) 

,  “  -  (aC^^  +  AC^^ ) 

L 

(38) 

(39) 

=  Xs+{Xs  -X4) 

(40) 
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A  trapezoidal  distribution  for  the  local  negative  contribution  to  the  body  loads  can  now  be  defined 
as  follows: 


x^<x<  x^: 


X-X^ 


(41) 


X^  <X<  Xgl 


Ac„ 


(42) 


X^  <  x<  Xy'. 


(43) 


These  local  load  components  are  added  to  the  total  local  body  loads  within  the  specified 
streamwise  locations.  If  there  is  insufficient  afterbody  length  to  allow  the  inclusion  of  the  entire 
negative  interference  load,  the  remainder  is  carried  over  to  the  leeward  tail  surfaces. 


2.3.2  Liftinu  Surface  Load  Changes 

When  a  missile  is  rolled  into  the  X  (<b  =  45  deg)  position,  the  lower  (windward  plane) 
and  upper  (leeward  plane)  lifting  surfaces  will  experience  different  aerodynamic  environments. 
The  normal  force  computed  by  the  AP95  code  is  the  total  for  all  fins,  both  windward  and  leeward, 
and  it  must  be  divided  between  the  two  planes.  In  order  to  accomplish  this  division,  the  approach 
previously  used  in  Reference  2  is  employed.  It  is  assumed,  in  agreement  with  slender  body 
theory,  that  at  an  AOA  of  zero,  the  normal  force  is  evenly  distributed  between  the  windward  and 
leeward  planes.  As  the  AOA  is  increased,  the  load  on  the  windward  surfaces  is  increased  linearly 
up  to  90  percent  of  the  total  at  a  =  65  deg.  This  percentage  remains  constant  up  to  a  =  90  deg. 
On  the  other  hand,  the  load  on  the  leeward  plane  decreases  in  a  corresponding  fashion  so  that  the 
total  of  the  two  is  not  changed.  Thus,  we  have 


a  <  65°:  [c„iy)]ymDWARD  =  [c„(y)](^0.5+0.4^ 


— 


(44) 

(45) 
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a  >65°: 


[^n  ()')]}raVDWlRD  (^•^[*'n(y)] 

(46) 

(47) 

A  further  point  to  consider  is  that  the  normal  forces  which  are  computed  by  the  AP95  code  are 
normal  only  to  horizontal  hfting  surfaces  in  the  <J>  =  0  deg  roll  orientation.  In  the  <I)  =  45  deg 
case,  the  loads  given  by  Equations  (32)-(35)  must  be  divided  by  cos  O  on  the  windward  side  and 
by  sin  0  on  the  leeward  side  to  give  locally  normal  values. 


The  portion  of  the  wing-tail  interference,  ,  which  is  allocated  to  the  tail  surfaces 

must  also  be  divided  between  the  windward  and  leeward  planes.  Once  again  rel5/ing  on 
methodology  developed  in  Reference  2,  this  division  is  accomplished  using  the  wing-tail 
interference  factors  defined  during  the  slender  body  derivation  of  this  term.  The  mathematics  of 
this  process  are  quite  involved  and  wiU  not  be  repeated  here.  For  more  information,  the  reader 
may  consult  the  cited  reference.  For  the  present  purpose,  if  ii  is  the  interference  factor  associated 
with  the  windward  plane  fin  and  u  the  one  with  the  leeward  plane  fin,  then 


WINDWARD 


(48) 


and 


LEEWARD 


\c  1 

1  ^T(V) 

Kh  +(47 

(49) 


This  brings  us  to  the  second  physical  phenomenon  which  is  not  expUcitly  modeled  in  the 
approach  outlined  above.  In  the  O  =  45  deg  roll  position,  there  is  a  “shadowing”  effect  on  the 
leeward  plane  lifting  surfaces  for  high  Mach  number  flows.  A  “shadowed  region”  exists  near  the 
root  chord  where  Newtonian  theory  predicts  a  pressure  coefficient  of  zero.  As  a  result,  the 
predicted  values  of  loads  on  the  leeward  plane  lifting  surfaces  to  be  presented  later  will  be  high, 
but  due  to  the  complexities  involved  in  trying  to  model  this  effect,  no  attempt  has  been  made  to 
do  so  at  this  time.  The  higher  the  Mach  number,  the  more  accurate  the  Newtonian  prediction  of 
zero  pressure  coefficient  in  the  shadowed  region  becomes.  This  therefore  means  that  the  current 
leeward  plane  loads  near  the  body  will  become  less  accurate  as  Mach  number  increases. 

The  final  phenomenon  which  is  believed  to  affect  the  load  distributions  but  is  not  yet  fully 
modeled  involves  the  stronger  influence  of  body  vortices  on  the  leeward  plane  fins  as  their 
location  is  moved  farther  aft  on  the  body.  In  the  development  of  the  semi-empirical  approach  of 
the  AP95,  the  body  vortex  effects  were  inherently  included  in  the  body  carryover  load.  Hence, 
the  body  carryover  load  was  lower  for  the  tail  fins  than  it  should  have  been  to  offset  the  fact  that 
the  vortices  actually  cause  the  leeward  plane  tail  fins  to  lose  more  normal  force  the  farther  aft  they 
are  located.  To  account  for  this  fact,  negative  contributions  from  the  adjusted  body-lifting  surface 
interference  loads,  as  described  in  Section  2.3.1,  which  cannot  be  distributed  on  the  afterbody  due 
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to  lack  of  sufficient  length  will  be  carried  over  to  the  leeward  plane  tail  fins.  They  wiU  be 
distributed  along  the  span  in  the  same  manner  as  the  linear  and  nonhnear  lifting  surface  loads. 


2.4  LOADS,  SHEAR,  AND  BENDING  MOMENTS 


The  primary  purpose  of  computing  load  distributions  over  the  components  of  a  missile 
body  is  to  allow  the  determination  of  shear  and  bending  moment  within  the  structural  members. 
The  loads  which  have  been  determined  so  far  are  non-dimensional  in  the  sense  that  they  represent 
a  normal  force  coefficient  per  unit  length.  Multiplying  them  by  the  dynamic  pressure  and  the 
reference  area  will  give  local  loads  in  pounds  per  foot.  These  loads  can  be  integrated  to  get  the 
local  shear  in  pounds  and  the  shear  can  in  turn  be  integrated  to  get  the  local  moment  in  pound- 
feet.  For  the  body,  utilizing  the  results  of  Equations  (2),  (6),  (7),  (8),  (17),  (18),  (19),  and  (30), 
we  have  for  the  dimensionalized  local  loads  w(x) 


w(x)„  =  (c„  )  +fc„  )  +Ac„  +Ac„  +(c„  ) 

^  \\  ”b  /L  \  JnL,ADJ  ^B(W)  ^b(T)  \  ”b  /T(y) 


QA 


REF 


where 


and 


Q  =  ^P^Ml 


Aref  ~  ^  ^ 


(50) 

(51) 

(52) 


We  can  then  integrate  to  get  the  shear,  V(x)b,  and  the  bending  moment,  BM(x)b. 


X 


0 

(53) 

X 

BM{x)^=jv{x)^dx 

(54) 

0 


The  procedure  for  the  lifting  surfaces  is  the  same  except  for  the  spanwise  rather  than  the 
streamwise  variation  of  the  loads.  For  the  wing,  using  the  nondimensional  load  from 
Equation  (22),  we  get 


w(y)w  =  [{c„(y))yv]QA  REF  •  (55) 

For  the  tail,  there  may  be  a  contribution  from  the  wing-tail  interference  given  by  Equation  (30)  in 
addition  to  the  nondimensional  load  from  Equation  (23). 

Hy)T  =  [(^"  (>'))7’  +  iy))nv)  ]Q^i^eF  (56) 
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Equations  (55)  and  (56),  as  presented  above,  are  valid  for  the  O  =  0  deg  roll  position.  For  the 
<I)  =  45  deg  case,  there  will  be  separate  loads  for  the  windward  and  leeward  sets  of  lifting 
surfaces.  These  separate  loads  can  be  computed  simply  by  substituting  the  windward  and 
leeward  values  of  the  (Cn(y))w,  (Cn(y))T,  and  (Cn(y))T(V)  contributions  given  in  Equations  (32)-(35). 
The  integration  for  shear  and  bending  moment  wiU  be  the  same  for  both  lifting  surfaces  and  for 
both  roll  orientations.  With  the  W/T  subscript  indicating  the  appropriate  value  for  either  the  wing 
or  tail,  we  have 


y 

0 

y 

^^{y)wiT  ~  • 

0 


(57) 

(58) 


All  of  the  above  integrations  are  performed  in  the  AP95  code  using  the  trapezoidal  rule. 


3.0  RESULTS  AND  DISCUSSION 


For  purposes  of  validation,  the  results  computed  by  the  new  Aeroprediction  methodology 
are  compared  to  Navier-Stokes*'’  computer  calculations  performed  by  Clint  Housh  at  the  Naval 
Air  Warfare  Center  Weapons  Division  (NAWCWPNS).  Distributions  of  aerodynamic  loads, 
shear,  and  bending  moment  were  generated  by  both  methods  for  the  SEASPARROW  air-to-air 
missile  at  Mach  numbers  of  1.50,  2.87,  and  4.60  and  angles  of  attack  of  10,  25,  and  40  deg.  The 
SEASPARROW,  shown  in  Figure  6,  with  its  two  sets  of  lifting  surfaces,  provides  an  opportunity 
to  exercise  aU  of  the  new  additions  to  the  code.  The  results  of  these  computations  are  presented 
and  discussed  below.  The  general  format  followed  is  to  handle  each  component  (body,  wing,  tail) 
separately.  First,  the  local  loads  are  plotted  on  three  separate  graphs  corresponding  to  the  three 
angles  of  attack.  AU  three  Mach  numbers  are  shown  on  the  same  plot.  The  shear  and  bending 
moment  are  then  presented,  in  turn,  in  groups  of  three  graphs  corresponding  to  the  angles  of 
attack.  This  sequence  is  then  repeated  for  the  next  missile  component.  In  aU  cases,  the  CFD 
results  are  plotted  together  with  the  Aeroprediction  results  for  comparison  purposes.  Before 
discussing  the  comparison  of  the  semiempMcal  and  CFD  results,  a  brief  discussion  of  the  CFD 
methodology  is  in  order. 


3.1  CFD  METHODOLOGY 

For  the  purposes  of  validation,  the  results  computed  by  the  new  Aeroprediction 
methodology  are  compared  to  Navier-Stokes^’  computer  calculations  performed  by  Clint  Housh 
at  the  Naval  Air  Warfare  Center  Weapons  Division  (NAWCWPNS).  The  configuration  chosen  to 
provide  the  validation  data  for  the  new  semiempirical  model  was  the  SEASPARROW  missile 
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configuration.  This  configuration  was  chosen  because  of  the  availability  of  experimental  force  and 
moment  data,”  which  was  used  to  provide  a  level  of  confidence  in  the  CFD  results.  The 
configuration  consisted  of  an  ogive-cylinder  body  with  a  fineness  ratio  of  18,  cruciform  wings, 
and  in-line  tails.  The  wings  had  a  trapezoidal  planform  with  a  leading-edge  sweep  of  45  deg  and  a 
diamond  airfoil  section.  The  tails  had  a  delta  planform  with  a  leading-edge  sweep  of  57  deg  and  a 
modified  diamond  airfoil  section.  Additional  geometrical  detail  can  be  found  in  Reference  11. 

Computational  grids  were  developed  for  the  SEASPARROW  geometry  in  both  the  “+” 
(roU  angle  0  =  0  deg)  configuration,  and  the  “x”  (roU  angle  O  =  45  deg)  configuration.  The  CFD 
grid-models  for  these  analyses  were  developed  using  the  Chimera  overset  mesh  approach.*^  The 
Chimera  overset  mesh  technique  is  a  domain  decomposition  procedure  where  a  configuration  is 
meshed  using  a  collection  of  overset  grids.  Grids  for  the  individual  components  of  a  configuration 
are  developed  independently  of  one  another  and  are  overset  upon  the  others.  Because  each 
component  grid  is  generated  independently  (i.e.,  body  grid,  wing  grid,  tail  grid),  portions  of  one 
grid  may  he  inside  the  physical  surfaces  of  another.  These  points  are  identified  and  excluded  from 
the  computations  using  the  Pegasus  computer  program.”  The  Pegasus  code  defines  regions 
where  holes  are  cut  by  one  grid  into  another  and  provides  the  interpolation  stencils,  used  by  the 
flow  solver,  at  the  various  grid  interfaces.  The  grids  of  the  individual  components  were 
developed  using  the  Gridgen*'*  computer  program. 

Viscous  grid  spacing  was  employed  over  all  of  the  sohd  surfaces  in  order  to  obtain 
accurate  solutions  to  the  Navier-Stokes  equations.  The  grid  spacing  normal  to  sohd  walls  was  set 
to  approximately  0.00004  body  diameters,  corresponding  to  y'*'  values  on  the  order  of  one.  The 
grid  model  for  the  “+”  configuration  was  made  up  of  7  separate  grids  totaling  2,131,950  grid 
points.  The  grid  model  for  the  “x”  configuration  has  5  grids  consisting  of  2,223,500  grid  points. 

The  OVERFLOW^”  computer  program  was  used  to  provide  solutions  to  the  thin-layer 
formulation  of  the  Reynolds-averaged  Navier-Stokes  equations.  Among  the  many  options 
available  in  the  code,  the  ARC3D-based”  imphcit  3-factor  diagonal  scheme  was  chosen  to 
advance  the  solution  in  time.  Conventional  central  differences  were  used  to  evaluate  the  inviscid 
and  viscous  terms  in  each  computational  direction.  Second  and  fourth  order  smoothing  were 
apphed  to  both  the  right  and  left  hand  sides.  Also,  for  the  energy  equation,  the  quantity  ph  was 
smoothed  rather  than  pe.  The  Baldwin-Barth  one-equation  turbulence  model‘s  was  used  for 
closure. 

The  boundary  condition  used  for  solid  surfaces  was  that  of  an  isothermal  no-slip  wall  with 
pressure  extrapolation.  The  downstream  overflow  boundary  was  set  using  simple  extrapolation. 
The  collapsed  face  from  the  tip  of  the  nose  to  the  upstream  boundary  plane  was  handled  by  setting 
the  boundary  value  to  an  average  of  all  surrounding  cells  by  varying  the  circumferential  index  and 
holding  the  axial  and  normal  indicies  constant.  On  the  symmetry  plane,  the  velocity  normal  to  the 
plane  was  set  to  zero.  The  upstream  boundary  and  far-field  boundary  values  were  set  using  a 
characteristic  condition  based  on  Riemann  invariants.  The  entire  flow  field  was  initialized  to  free 
stream  values  appropriate  to  the  experimental  flow  conditions  described  in  Reference  1 1. 
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The  solutions  were  run  on  a  Silicon  Graphics  Incorporated  (SGI)  Power  Challenge  Array 
using  a  single  R8000  processor.  The  OVERFLOW  program  was  run  with  the  namelist  option: 
INCORE=FALSE.  In  this  mode,  memory  requirements  are  set  by  the  largest  individual  grid. 
Grid,  flowfield  and  interpolation  arrays  are  read  from  scratch  files,  one  (or  more)  iterations  are 
taken,  and  the  flowfield  is  written  back  out  to  disk.  The  memory  required  to  run  these  cases  on 
the  SGI  machine  was  on  the  order  of  9.4  Mwords  (8-byte  words).  The  processing  time  used  by 
the  computer,  CPU  time,  was  approximately  62.33  nsec/grid  point/iteration.  A  typical  solution 
took  3000  iterations  to  converge  the  viscous  and  inviscid  forces.  Local  time  stepping  was  used  to 
accelerate  the  solutions  to  steady  state. 


Solutions  have  been  obtained  for  the  SEASPARROW  missile  configuration  at  three 
supersonic  Mach  numbers  (M  =  1.5,  2.87  and  4.63),  two  roll  angles  (<h  =  0  and  45  deg),  and 
three  different  angles  of  attack  (a  =  10,  25,  and  40  deg).  Table  1  shows  a  summary  of  the 
solutions  obtained  and  the  percent  difference  between  experimental  and  Navier-Stokes  computed 
lift  and  drag  coefficients.  Also  shown  is  the  difference  in  center  of  pressure  location  between 
experimental  and  computed  results  expressed  as  a  percentage  of  the  overall  length  of  the  missile. 
The  comparisons  between  the  experimental  data  and  the  computed  results  allow  us  to  evaluate  the 
accuracy  of  the  solution  and,  to  some  extent,  validate  the  solution  methodology. 

TABLE  1.  SPARROW  HI  RUN  MATRIX 


Mach 

Number 

AOA 

a 

Roll 

O 

CL 
%6iff 
exp  vs.  CFD 

CD 
%diff 
exp  vs.  CFD 

Xcpdiff 
%  body  length 
exp  vs.  CFD 

10 

0 

1.8 

6.1 

0.78 

25 

0 

1.1 

2.6 

0.57 

1.50 

40 

0 

-0.7 

0.4 

0.07 

2.87 

10 

0 

2.1 

5.9 

1.44 

2.87 

25 

0 

-1.5 

-0.7 

0.07 

2.87 

40 

0 

-1.4 

-2.9 

0.38 

10 

0 

0 

nmgiiiiiiiiii 

1.13 

25 

0 

1.2 

0.35 

1  4.60 

40 

0 

-6.1 

2.42 

■11^  nH 

10 

45 

2.4 

■IQQiilll 

0.85 

25 

45 

0.4 

0.44 

1.50 

40 

45 

-1.5 

-0.09 

2.87 

10 

-1.6 

MBMM 

2.87 

25 

-3.0 

2.87 

40 

45 

-0.9 

0.56 

4.60 

10 

45 

-5.6 

6.7 

1.58 

4.60 

25 

45 

-0.6 

3.6 

1.21 

4.60 

40 

45 

0.5 

1.8 

0.82 

As  shown  in  the  table,  the  average  percent  difference  between  experimental  and  CFD 
values  of  lift  coefficient,  for  all  of  the  solutions  obtained,  was  shghtly  less  than  2  percent.  This  is 
good  agreement  between  the  experimental  and  CFD  results.  A  couple  of  the  solutions  did  not 
predict  the  lift  coefficient  as  well:  the  M  =  4.6,  a  =  40  deg,  <!>  =  0  deg  and  the  M  =  4.6, 
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a  =  10  deg,  <I)  =  45  deg  cases.  In  these  cases,  the  percent  difference  of  lift  coefficient  was  on  the 
order  of  6  percent.  For  these  two  cases  the  poor  prediction  of  the  magnitude  of  lift  was  also 
apparent  in  the  center  of  pressure  prediction. 

Also  shown  in  the  table,  the  computed  center  of  pressure  was  shghtly  aft  of  the 
experimental  center  of  pressure  for  ah  of  the  cases  run  (except  the  M  =  1.5,  a  =  40  deg, 
0  =  0  deg  case).  The  average  difference  between  the  experimental  and  computed  center  of 
pressure  results  was  on  the  order  of  1  percent  of  the  body  length  for  the  solutions  run.  In 
addition,  at  each  Mach  number,  the  percent  difference  of  drag  coefficient  between  experimental 
and  computed  results  was  greater  at  the  low  angles  of  attack.  This  is  not  surprising  considering 
the  large  increase  in  drag  coefficient  with  angle  of  attack.  If,  for  example,  the  experimental  and 
computed  drag  coefficients  are  separated  by  a  constant  6  regardless  of  angle  of  attack,  that  8  is  a 
larger  percentage  of  a  smaher  number. 

Ah  in  ah,  the  CFD  results  compared  favorably  with  the  experimental  data.  The  good 
agreement  over  the  wide  angle  of  attack  range  shows  the  solution  methodology  to  be  accurate  in 
predicting  both  linear  and  nonlinear  aerodynamics.  The  close  agreement  in  the  center  of  pressure 
prediction  show  the  capabhity  of  the  solution  methodology  to  accurately  model  the  load 
distribution.  The  good  agreement  validates  the  solution  methodology  and  confirms  the  accuracy 
of  the  solutions  obtained. 


3.2  RESULTS  FOR  <D  =  0  DEG  ROLL  POSITION 

Comparisons  of  the  body  results  from  the  new  Aeroprediction  methodology  and  the  CFD 
computations  are  shown  in  Figures  7  through  9.  Figure  7  shows  the  body  loads.  Agreement  is,  in 
general,  quite  good.  However,  even  after  the  adjustments  made  to  the  nose,  there  is  stUl  some 
slight  overprediction  of  the  loads  in  this  region.  In  addition,  for  the  M  =  4.60,  a  =  40  deg  case, 
there  is  some  “smearing  out”  of  the  Aeroprediction  body  interference  loads  in  the  vicinity  of  the 
wing.  This  could  be  the  result  of  shock  interaction  effects  which  are  not  modeled  in  the  code. 
Figure  8  shows  the  comparisons  for  shear  force  at  the  same  corresponding  angles  of  attack.  The 
Aeroprediction  code  does  a  very  good  job  with  some  overprediction  for  the  a  =  10  deg  cases 
relative  to  the  CFD  results.  Note  that  the  shear  results  are  very  good  for  the  M  =  4.60, 
a  =  40  deg  case  cited  above.  The  under  and  over  prediction  regions  for  local  loads  tend  to  cancel 
each  other  out  when  they  are  integrated  to  get  the  shear.  The  bending  moment  results  for  the 
body  are  plotted  in  Figure  9.  Once  again,  agreement  is  quite  good,  with  some  slight 
overprediction  at  low  a. 

Figures  10-12  present  the  results  for  the  wing.  The  local  loads  are  shown  in  Figure  10. 
Agreement  with  the  CFD  results  is  excellent  in  most  cases  with  the  only  exceptions  being  at 
M  =  4.60  near  the  body  where  the  loads  are  overpredicted.  It  is  possible  that  there  are  some 
shock  interaction  effects  which  are  not  accounted  for  in  this  region.  Wing  shear  results  are 
plotted  in  Figure  1 1  and  they  are  in  very  good  agreement  with  the  CFD  results  with  only  small 
overpredictions  at  a  =  10  deg.  Bending  moments,  shown  in  Figure  12,  are  also  predicted  quite 
well  with  some  slight  underprediction  at  high  a. 
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The  tail  comparisons  are  shown  in  Figures  13-15.  The  loads,  presented  in  Figure  13, 
while  not  in  quite  as  close  agreement  as  for  the  wing,  are  still  good,  with  the  region  near  the  wing 
root  at  high  Mach  number  once  again  being  overpredicted.  Figure  14  shows  the  corresponding 
shear  loads.  The  correlation  with  the  CFD  results  is  once  again  quite  good  with  a  small  degree  of 
overprediction,  especially  at  a  =  40  deg.  The  bending  moment  results,  shown  in  Figure  15,  are  in 
very  good  agreement. 


3.3  RESULTS  FOR  0  =  45  DEG  ROLL  POSITION 

Figure  16  shows  the  comparisons  for  the  body  loads.  At  M  =  1.50,  the  Aeroprediction 
results  exhibit  good  agreement  with  the  CFD  computations.  However,  at  M  =  2.87  and  to  an 
even  greater  extent  at  M  =  4.60,  the  CFD  predictions  show  a  sharp  rise  for  the  body  loads  in  the 
vicinity  of  the  wings  followed  by  a  rapid  decrease  below  body  alone  values  behind  the  wings. 
This  behavior  is  only  partially  captured  in  the  Aeroprediction  results  which  show  a  smoother  rise 
with  a  lower  peak  and  a  lesser  decline  behind  the  wings.  As  discussed  in  Section  2.1.3, 
modifications  were  made  to  the  body  carryover  methodology  for  the  $  =  45  deg  roll  position  to 
improve  the  modeling  of  this  phenomenon.  These  changes  have  resulted  in  a  great  improvement 
in  the  code’s  performance  in  this  region,  but  the  sharp  surface  pressure  variations  are  still  not 
completely  modeled.  The  physical  mechanism  which  gives  rise  to  the  pressure  spike  and  the 
following  overexpansion  is  not  completely  understood  at  the  present  time  and,  since  it  is  most 
likely  associated  with  strong  shock  interactions  on  the  surface,  it  is  considered  beyond  the  scope 
of  the  present  work  to  try  and  develop  a  more  rigorous  model.  Figure  17  shows  the  results  for 
the  shear  force  distributions.  The  agreement  with  the  CFD  results  is  quite  good,  especially 
considering  the  body  load  variations.  The  body  bending  moments,  shown  in  Figure  18,  exhibit 
even  better  agreement. 

The  spanwise  load  distributions  for  the  windward  plane  wings  are  plotted  in  Figure  19. 
The  agreement  between  the  Aeroprediction  and  CFD  results  is  very  good  with  only  some  minor 
deviations  for  M  =  4.60.  The  shear  loads  are  shown  in  Figure  20.  Again,  the  agreement  is  quite 
close  with  the  Aeroprediction  results  being  slightly  high  for  M  =  4.60  at  the  two  lower  angles  of 
attack.  The  bending  moments,  presented  in  Figure  21,  show  a  general  tendency  to  be  slightly  low 
relative  to  the  CFD  resists  but  overall  agreement  is  good. 

The  load  distributions  for  the  leeward  wings  are  shown  in  Figure  22.  It  can  be  seen  that 
the  “shadowing”  effect  discussed  above  is  quite  evident  for  the  M  =  2.87  and  M  =  4.60  cases. 
The  CFD  results  predict  a  large  pressure  loss  on  the  wings  in  the  shadowed  region  near  the  body 
which  is  not  modeled  in  the  Aeroprediction  code.  Otherwise,  agreement  is  very  good.  The  shear 
distributions  are  shown  in  Figure  23,  and,  as  might  be  expected,  there  is  significant  overprediction 
at  the  higher  Mach  numbers  relative  to  the  CFD  results,  but  the  results  are  still  surprisingly  good. 
The  same  comments  apply  to  the  bending  moments  which  are  presented  in  Figure  24. 

Figure  25  shows  the  windward  tail  loads.  The  agreement  between  the  two  approaches  is 
quite  good  except  for  the  M  =  4.60  case  at  a  =  25  deg.  The  CFD  results  show  a  loss  of  pressure 
on  the  inboard  section  of  the  tail  surfaces  which  may  be  associated  with  a  shock  or  vortex 
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impingement  which  is  not  modeled  by  the  Aeroprediction  code.  With  this  exception,  the  shear 
and  bending  moment  distributions,  shown  in  Figures  26  and  27,  respectively,  show  very  good 
agreement. 

The  leeward  tail  loads  are  plotted  in  Figure  28.  As  for  the  leeward  wings,  the 
“shadowing”  effect  is  very  obvious  at  the  two  higher  Mach  numbers.  The  shear  loads,  shown  in 
Figure  29,  and  the  bending  moments,  shown  in  Figure  30,  tend  to  overpredict  relative  to  the  CFD 
results  because  of  this  effect.  It  should  be  noted  that  by  modifying  the  body  carryover  prediction 
methodology  for  the  O  =  45  deg  case  and  applying  any  excess  negative  contributions  to  the 
leeward  tail,  the  results  for  the  leeward  tail  have  been  improved  substantially  over  their  original 
values.  However,  because  of  the  complex  nature  of  leeside  flow  fields  at  high  angles  of  attack,  it 
is  extremely  difficult  to  develop  empirical  models  that  will  give  truly  accurate  results  in  these 
regions. 


3.4  RESULTS  FOR  WING  DEFLECTION 

As  a  final  check  on  the  new  additions  to  the  code,  a  case  was  run  at  M  =  2.87  and 
a  =  10  deg  for  a  wing  deflection  of  20  deg.  Only  the  ^)  =  0  deg  roll  position  was  considered. 
The  resulting  loads,  shear,  and  bending  moment  for  the  body  are  shown  in  Figure  31.  These  loads 
reflect  contributions  from  body- wing  interference  factors,  Kb(W)  and  kB(W),  associated,  respectively, 
with  the  total  missile  angle  of  attack  and  with  the  wing  deflection  angle.  As  should  be  expected, 
the  curves  for  the  imdeflected  and  deflected  case  are  identical  up  to  the  location  of  the  wing 
where  the  body  carryover  load  from  the  deflected  wing  is  much  higher  than  for  the  undeflected 
case.  There  is  also  a  slightly  lower  load  in  the  region  of  the  tail  when  the  wing  is  deflected 
because  of  the  more  negative  wing-tail  interference  effect. 

The  curves  for  the  loads,  shear,  and  bending  moment  of  the  wing  are  presented  in 
Figure  32.  The  results  are  once  again  as  would  be  expected  with  higher  values  for  aU  parameters 
when  the  wing  is  deflected. 

Finally,  Figure  33  show  the  loads,  shear,  and  bending  moment,  respectively,  for  the  tail. 
Note  that,  in  this  case,  because  of  the  more  negative  effect  of  the  wing-tail  interference,  the  loads 
on  the  tail  are  lower  when  the  wing  is  deflected.  The  lower  loads  translate  into  lower  values  for 
the  shear  and  bending  moment  as  well. 


4.0  SUMMARY 


New  methodology  has  been  added  to  the  NSWCDD  Aeroprediction  code  to  compute 
both  the  linear  and  nonlinear  contributions  to  surface  loads  and  to  distribute  these  loads  over  the 
body  and  lifting  surfaces.  The  new  methods  include  models  to  simulate  hfting  surface  interference 
effects  on  the  body  for  both  the  O  =  0  deg  and  3)  =  45  deg  roll  positions.  At  O  =  45  deg,  load 
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distributions  are  computed  for  both  the  windward  and  leeward  lifting  surfaces.  Wing-tail 
interference  loads  are  now  separated  into  separate  components  which  are  distributed  along  the 
body  in  the  vicinity  of  the  tail  as  well  as  across  the  tail  surfaces  themselves.  All  loads  are 
integrated  to  give  shear  and  bending  moment  distributions  which  may  be  used  by  structural 
engineers  in  the  design  process. 

As  an  aid  in  developing  the  new  techniques  and  for  validation  purposes,  results  from 
Navier-Stokes  CFD  computations  for  the  SEASP ARROW  air-to-air  missile  were  utilized.  In 
general,  the  Aeroprediction  methodology  proved  to  be  quite  successful  at  modeling  the 
distribution  of  surface  loads  when  compared  to  the  CFD  results.  The  only  exceptions  occurred 
for  the  <5  =  45  deg  roll  position  where  the  loads  on  leeside  lifting  surfaces  (especially  tail  surfaces) 
tended  to  be  overpredicted,  and  the  exaggerated  pressure  variations  on  the  body  in  the  vicinity  of 
lifting  surfaces  were  not  fully  captured. 

With  these  new  additions,  the  Aeroprediction  code  should  prove  to  be  a  much  more  useful 
tool  to  the  structural  engineer  who  is  interested  in  performing  preliminary  structural  analyses  of 
several  competing  missile  designs. 
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AOA 

AP95 

AP98 

CFD 

NAWCWPNS 

NSWCDD 

Ap 

Aref 

br 

bw 

BM(x)b 

BM(y)w/T 

c(y) 

c(y)T 

C(y)w 

)mAX 

)mIN 


6.0  SYMBOLS  AND  DEFINITIONS 

Angle  of  Attack 

Aeroprediction  1995 

Aeroprediction  1998 

Computational  Fluid  Dynamics 

Naval  Air  Warfare  Center,  Weapons  Division 

Naval  Surface  Warfare  Center,  Dahlgren  Division 

Planform  area  of  the  body  in  the  crossflow  plane  (ft^) 

Reference  area  (maximum  cross-sectional  area  of  body,  if  a  body  is 
present,  or  planform  area  of  wing  if  wing  alone)  (ft^) 

Tail  span  (not  including  body)  (ft) 

Wing  span  (not  including  body)  (ft) 

Body  bending  moment  as  a  function  of  stream  wise  coordinate  (Ib-ft) 
Wing  or  tail  bending  moment  as  a  function  of  spanwise  coordinate  (Ib-ft) 
Local  chord  of  lifting  surface  as  a  function  of  spanwise  coordinate  (ft) 
Local  chord  of  tail  as  a  function  of  spanwise  coordinate  (ft) 

Local  chord  of  wing  as  a  function  of  spanwise  coordinate  (ft) 

Local  normal  force  coefficient  of  the  body 

Linear  component  of  body  local  normal  force  coefficient 

Maximum  value  of  local  body  normal  force  coefficient  in  trapezoidal 
distribution  of  positive  hfting  surface  interference  contribution 

Minimum  value  of  local  body  normal  force  coefficient  in  trapezoidal 
distribution  of  negative  hfting  surface  interference  contribution 
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)nl 

Nonlinear  component  of  body  local  normal  force  coefficient 

V  ”fi  )nl,adj 

Nonlinear  component  of  body  local  normal  force  coefficient  adjusted  for 

decreased  nose  loads 

Component  of  body  local  nonlinear  normal  force  coefficient,  downstream 

of  nose,  added  to  make  up  for  decrease  of  normal  force  on  nose 

[(^”*  )wt]rfiAW 

Component  of  body  local  nonlinear  normal  force  coefficient, 

downstream  of  nose,  added  to  make  up  for  loss  of  normal  force  in 
transition  region 

)m:  ];c-£ 

Nonlinear  component  of  body  local  normal  force  coefficient  at  end  of 

nose 

(c  ) 

\  /T(y) 

Component  of  body  local  normal  force  coefficient  due  to  wing-tail 

interference 

Linear  component  of  local  normal  force  coefficient  on  lifting  surfaces  as  a 
function  of  spanwise  coordinate 

Linear  component  of  local  normal  force  coefficient  on  tail  as  a  function 

of  spanwise  coordinate 

Linear  component  of  local  normal  force  coefficient  on  wing  as  a  function 

of  spanwise  coordinate 

\pn  (y ^\ieEWARD 

Portion  of  total  local  normal  force  coefficient  assigned  to  the  leeward 
lifting  surfaces 

(c„(y))j. 

Local  normal  force  coefficient  on  tail  as  a  function  of  spanwise 

coordinate 

ipn  (y ))7’(7) 

Portion  of  local  normal  force  coefficient  on  lifting  surfaces  attributed  to 
wing-tail  interference 
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WlNDWAm) 


Cn 

.  ^  Jl 


)l 

)w 


Local  normal  force  coefficient  on  wing  as  a  function  of  spanwise 
coordinate 

Portion  of  total  local  normal  force  coefficient  assigned  to  the  windward 
lifting  surfaces 

Total  normal  force  coefficient 
Nonnal  force  coefficient  of  the  body 

Total  negative  contribution  of  lifting  surface  interference  to  body  normal 
force 

Body  normal  force  coefficient  contribution  from  region  A  of  trapezoidal 
distribution 

Body  normal  force  coefficient  contribution  fi’om  region  B  of  trapezoidal 
distribution 

Body  normal  force  coefficient  contribution  from  region  C  of  trapezoidal 
distribution 

Linear  component  of  streamwise  variation  of  body  normal  force 
coefficient  (ff^) 

Nonlinear  component  of  streamwise  variation  of  body  normal  force 
coefficient  (ff*) 

Linear  component  of  the  body  normal  force  coefficient 
Nonlinear  component  of  the  body  normal  force  coefficient 
Increment  to  normal  force  coefficient  on  body  due  to  presence  of  tail 
Increment  to  normal  force  coefficient  on  body  due  to  presence  of  wing 
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®=45’ 


T(B) 


Q. 


BODY 


LEEWARD 


TAIL 


iwflVDWARD 


Cp 

Q 

iCr)r 


{CX 

C., 

Db 


e 

ii.it 


Contribution  to  body  normal  force  due  to  presence  of  a  lifting  surface  at 
<1>  =  45  deg  roll  position 

Normal  force  coefficient  of  tail  in  presence  of  body 

Component  of  normal  force  coefficient  in  tail  region  due  to  interference 
from  wing 

Portion  of  total  wing-tail  interference  term  assigned  to  body 

Portion  of  tail  component  of  wing-tail  interference  term  assigned  to  the 
leeward  tail  surfaces 

Portion  of  total  wing-tail  interference  term  assigned  to  tail 

Portion  of  tail  component  of  wing-tail  interference  term  assigned  to  the 
windward  tail  surfaces 

Normal  force  coefficient  of  the  wing  in  the  presence  of  the  body 

Change  of  wing  normal  force  coefficient  with  angle  of  attack  (rad'*) 

Pressure  coefficient 

Root  chord  of  lifting  surface  (ft) 

Root  chord  of  tail  (ft) 

Root  chord  of  wing  (ft) 

Crossflow  drag  coefficient 

Diameter  of  body  (ft) 
internal  energy  (ft^/sec^) 

Interference  factors  in  wing-tail  interference  calculations  for  windward 
and  leeward  plane  wings,  respectively 
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h 

k 

kB(W) 

kw(B) 

K 


Kb(W) 

AKb(w) 

\Kb(W)  ] 

MJN 

1^=0" 

jitf/iV, 

Kw(b) 


Lnose 

M 

Ml 

Moo 

Poo 

Q 

r 

Tvort 


local  enthalpy  (ftVsec^) 

Ratio  of  normal  force  contribution  of  a  deflected  lifting  surface  in  the 
presence  of  a  body  to  that  of  the  lifting  surface  alone  at  a  =  0  deg 

Ratio  of  additional  body  normal  force  coefficient  due  to  the  presence  of  a 
deflected  wing  to  that  of  the  wing  alone  at  a  =  0  deg 

Ratio  of  normal  force  contribution  of  a  deflected  wing  in  the  presence  of 
a  body  to  that  of  the  wing  alone  at  a  =  0  deg 

Ratio  of  lifting  surface  normal  force  coefficient  in  the  presence  of  a  body 
to  that  of  the  lifting  surface  alone  at  6  =  0  deg 

Ratio  of  additional  body  normal  force  coefficient  in  the  presence  of  a 
wing  to  that  of  the  wing  alone  at  5  =  0  deg 

Difference  between  minimum  Kb(W)  value  at  O  =  0  deg  and  d>  =  45  deg 
Minimum  value  of  Kb(W) 

Minimum  value  of  Kb(w)  at  the  <!>  =  0  deg  roll  position 

Minimum  value  of  Kb(W)  at  the  C>  =  45  deg  roll  position 

Ratio  of  wing  normal  force  coefficient  in  the  presence  of  a  body  to  that 
of  the  wing  alone  at  6  =  0  deg 

Length  of  the  nose 

Mach  number 

Local  Mach  number 

Freestream  Mach  number 

Freestream  pressure  (Ib/ft^) 

Dynamic  pressure  (Ib/ft^) 

Local  body  radius  (ft) 

Vortex  radius  (ft) 
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VC(y) 

V(x)b 

V(y)w/r 

w(x)b 

w(y)T 

w(y)w 

X 

XlE 

XtE 

X2,  X3,  X4 

X5,  X6,  X? 

y 

ybBT 

yvORT 

a 

p 

Y 

5 

11 


Vortex  strength-local  chord  parameter 

Body  shear  as  a  function  of  stream  wise  coordinate  (lb) 

Wing  or  tail  shear  as  a  function  of  spanwise  coordinate  (lb) 

Local  body  load  as  a  function  of  streamwise  coordinate  (Ib/ft) 

Local  tail  load  as  a  function  of  spanwise  coordinate  Gb/ft) 

Local  wing  load  as  a  function  of  spanwise  coordinate  Gb/ft) 

Streamwise  coordinate  (ft) 

Streamwise  location  of  leading  edge  (ft) 

Streamwise  location  of  trailing  edge  (ft) 

Streamwise  coordinates  used  to  define  the  trapezoidal  distribution  of  the 
positive  body-hfting  surface  interference  region  (ft) 

Streamwise  coordinates  used  to  define  the  trapezoidal  distribution  of  the 
negative  body-lifting  surface  interference  region  (ft) 

Spanwise  lifting  surface  coordinate  (ft) 

Spanwise  distance  from  path  of  vortex  shed  from  wing  to  local  point  on 
tail  (ft) 

Spanwise  distance  from  centerline  of  body  to  path  of  vortex  shed  from 
wing 

Angle  of  attack  (rad,  deg) 

Specific  heat  ratio  of  air 

Deflection  angle  of  control  surface  (rad,  deg) 

Parameter  used  in  viscous  crossflow  theory  for  nonlinear  body  normal 
force  (in  this  context,  it  is  the  ratio  of  the  normal  force  of  a  circular 
cylinder  of  given  length-to-diameter  ratio  to  that  of  a  cylinder  of  infinite 
length) 
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p  density  (slugs/ft^) 

O  Roll  position  of  missile  (<J>  =  0  deg  corresponds  to  fins  in  the  plus  (+) 

orientation) 
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HGURE  4.  GEOMETRY  OF  WING-TAIL  INTERFERENCE 
DISTRIBUTION  MODEL 
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M=4.60,AP 


HGURE  7.  BODY  LOADING  AT  4)  =  0°:  (a)  ALPHA  =  10° 
(b)  ALPHA  =  25°,  AND  (c)  ALPHA  =  40° 
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x/L 


(a)  a  =10“ 


x/L 


(b)  a  =  25“ 


(c)  a  =  40“ 


FIGURE  8.  BODY  SHEAR  AT  O  =  0“:  (a)  ALPHA  =  10' 
(b)  ALPHA  =  25“,  AND  (c)  ALPHA  =  40° 
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•M=4.60,AP 


FIGURE  10.  WING  LOADING  AT  «>  =  0°:  (a)  ALPHA  =  10°, 
(b)  ALPHA  =  25°,  AND  (c)  ALPHA  =  40° 
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Y/(B/2) 


(b)  a  =  25° 


0  0.2  0.4  0.6  0.8  1 

Y/(B/2) 

(c)  a  =  40° 


FIGURE  1 1 .  WING  SHEAR  AT  O  =  0°:  (a)  ALPHA  =  10 
(b)  ALPHA  =  25°,  AND  (c)  ALPHA  =  40° 
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Y/(B/2) 

(a)  a  =  10° 


0  0.2  0.4  0.6  0.8  1 

Y/(B/2) 

(b)a  =  25° 


0  0.2  0.4  0.6  0.8  1 

Y/(B/2) 

(c)a  =  40° 


HGURE  12.  WING  BENDING  MOMENT  AT  O  =  0°:  (a)  ALPHA  =  10°, 
(b)  ALPHA  =  25°,  AND  (c)  ALPHA  =  40° 
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'M=4.60,AP 


HGURE  13.  TAIL  LOADING  AT  O  =  0°:  (a)  ALPHA  =  10°, 
(b)  ALPHA  =  25°,  AND  (c)  ALPHA  =  40° 
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BENDING  MOMENT  (LB-FT)  BENDING  MOMENT  (LB-FO  BENDING  MOMENT  (LB*FT) 
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Y/(B/2) 

(c)  a  =  40° 

HGURE  15.  TAIL  BENDING  MOMENT  AT  $  =  0°:  (a)  ALPHA  =  10°, 
(b)  ALPHA  =  25°,  AND  (c)  ALPHA  =  40° 
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FIGURE  16.  BODY  LOADING  AT  O  =  45°:  (a)  ALPHA  =  10' 
(b)  ALPHA  =  25°,  AND  (c)  ALPHA  =  40 
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HGURE  17.  BODY  SHEAR  AT  O  =  45°:  (a)  ALPHA  =  10°, 
(b)  ALPHA  =  25°,  AND  (c)  ALPHA  =  40° 
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- M=4.60,AP 

■  M=4.60,CFD 


FIGURE  19.  WINDWARD  WING  LOADING  AT  4)  =  45°:  (a)  ALPHA  =  10°, 
(b)  ALPHA  =  25°,  AND  (c)  ALPHA  =  40° 
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M=4.60,AP 


(a)  a  =10“ 


HGURE  21.  WINDWARD  WING  BENDING  MOMENT  AT  ®  =  45°:  (a)  ALPHA  =  10°, 
(b)  ALPHA  =  25°,  AND  (c)  ALPHA  =  40° 
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■M=4.60,AP 


HGURE  22.  LEEWARD  WING  LOADING  AT  O  =  45°:  (a)  ALPHA  =  10°, 
(b)  ALPHA  =  25°,  AND  (c)  ALPHA  =  40° 
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HGURE  23.  LEEWARD  WING  SHEAR  AT  <0  =  45“:  (a)  ALPHA  =  10°, 
(b)  ALPHA  =  25°,  AND  (c)  ALPHA  =  40° 
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-M=4.60,AP 

M=4.60,CFD 


Y/(B/2) 
(a)  a  =  10° 


HGURE  24.  LEEWARD  WING  BENDING  MOMENT  AT  O  =  45°:  (a)  ALPHA  =  10°, 
(b)  ALPHA  =  25°,  AND  (c)  ALPHA  =  40° 
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(c)  a  =  40° 


FIGURE  26.  WINDWARD  TAIL  SHEAR  AT  <E>  =  45°:  (a)  ALPHA  =  10° 
(b)  ALPHA  =  25°,  AND  (c)  ALPHA  =  40° 
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HGURE  27.  WINDWARD  TAIL  BENDING  MOMENT  AT  O  =  45°:  (a)  ALPHA  =  10°, 
(b)  ALPHA  =  25°.  AND  (c)  ALPHA  =  40° 
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M=4.60,AP 


(a)  a  =10° 


(c)  a  =  40° 


HGURE  28.  LEEWARD  TAIL  LOADING  AT  «»  =  45°:  (a)  ALPHA  =  10°, 
(b)  ALPHA  =  25°,  AND  (c)  ALPHA  =  40° 
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FIGURE  30.  LEEWARD  TAIL  BENDING  MOMENT  AT  a>  =  45°:  (a)  ALPHA  =  10°, 
(b)  ALPHA  =  25°,  AND  (c)  ALPHA  =  40° 
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M=2.87,8=0'’ 


nGURE31.  EFFECT  OF  WING  DEFLECTION  AT 
a  =  10°,  0)  =  0°:  (a)  BODY  LOADING,  (b)  BODY  SHEAR, 
AND  (c)  BODY  BENDING  MOMENT 
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M=2.87,6  =  0'> 


FIGURE  32.  EFFECT  OF  WING  DEFLECTION  AT 
a  =  10°,  =  0°:  (a)  WING  LOADING,  (b)  WING  SHEAR, 

AND  (c)  WING  BENDING  MOMENT 
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- M=2.87,6  =  0“ 

■  M=2.87,8  =  20° 


Y/(B/2) 


(c)  TAIL  BENDING  MOMENT 


nGURE33.  EFFECT  OF  WING  DEFLECTION  AT 
a  =  10°,  $  =  0°:  (a)  TAIL  LOADING,  (b)  TAIL  SHEAR, 
AND  (c)  TAIL  BENDING  MOMENT 
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